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Abstract 

Two applications of the Noether method for fluids and plasmas are presented based on the Euler- 
Lagrange and Euler-Poincare variational principles, which depend on whether the dynamical fields 
are to be varied independently or not, respectively. The relativistic cold laser-plasma equations, 
describing the interaction between an intense laser field with a cold relativistic electron plasma, 
provide a useful set of equations amenable to both variational formulations. The derivation of con- 
servation laws by Noether method proceeds from the Noether equation, whose form depends on the 
variational formulation used. As expected, the expressions for the energy-momentum conservation 
laws are identical in both variational formulations. The connection between the two Lagrangian 
densities is shown to involve the mass conservation and Lin constraints associated with the cold 
relativistic electron fluid. 

PACS numbers: 03.50.-z, 52.35.Mw 
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I. INTRODUCTION 

The Noether method is used in many areas of physics to relate symmetry properties 
of Lagrangians to important conservation laws for dissipationless dynamical equations such 
as energy, momentum, and wave action. Most conventional applications of the Noether 
method in classical and quantum field theories focus on a variational formulation that 

considers independent variations of all dynamical fields. The purpose of this paper is to show 
how the Noether method can be applied to derive the energy-momentum conservation law 
for a set of dissipationless plasma-fluid equations from two different variational formulations, 
which differ by their treatment of field variations. 

We begin with an expression for a typical action functional in field theory Q]: 



where the multi-component field ^> a (x, t) represents the state of the system at position x and 
time t (the index a is used to represent a field component) and we assume that the multi- 
component field is constrained to the space-time domain T> in the course of its evolution. 
In the present Section, a covariant notation is used in terms of the Minkowski metric tensor 
_ diag(l, — 1, — 1, — 1), the infinitesimal space-time volume is d A x = cdtd 3 x, and <9 M = 
(c~ l d t , V) is the covariant space-time gradient. The field equation for the dissipationless 
time evolution of the multi-component field V( x > t) is derived from a variational principle 
based on variations of the action functional 



where the variations 5tp a are assumed to vanish on the boundary dT> of the space-time 
domain T>. 

The stationarity of the action leads to the variational principle 



The outcome of the variational principle Q is a set of partial differential equations involving 
the Lagrangian density £ and its partial derivatives while the specific form of these equations 
depends on whether the variations 5ip a are independent or not. These two cases lead to the 
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Euler-Lagrange equations |3j and the Euler-Poincare equations [5|, respectively. As a result 
of the variational principle Q, the variation 5C of the Lagrangian density becomes the 
Noether equation 

8C = d,A», (4) 

where the four-density A M is expressed in terms of partial derivatives of the Lagrangian 
density C. The Noether method associates symmetries of the Lagrangian density C with 
local conservation laws as will be shown below. 

The remainder of the paper is organized as follows. In Sec. II, the standard application 
of the Noether method based on an Euler-Lagrange variational principle is reviewed. A 
general expression for the Euler-Lagrange form of the energy-momentum conservation law is 
derived by Noether method based on symmetries of the Euler-Lagrange Lagrangian density 
with respect to space-time translations. Some of the invariance properties of the energy- 
momentum conservation law are also investigated within the context of the symmetry of the 
energy-momentum tensor. 

In Sec. Ill, the relativistic cold laser-plasma equations p describing the interaction be- 
tween an intense laser field with a cold relativistic electron plasma in the presence of a neu- 
tralizing fixed-ion background. This example provides a useful set of equations amenable 
to both Euler-Lagrange and Euler-Poincare variational formulations and each application of 
the Noether method leads to the same energy-momentum conservation law. In the Euler- 
Lagrange variational formulation, the dynamical fields (ip, A; n, ip, a, j3) include the elec- 
tromagnetic potentials (<p, A), the electron density n, and the Clebsch potentials (ip, a, (3) 
used to represent the relativistic electron canonical momentum P = Vip + a V/3 (in the 
present work, the curl-free case P = discussed in Ref. jfj is generalized by considering 
V X P = Va X V/3 7^ 0). In the Euler-Poincare variational formulation, on the other hand, 
the Clebsch potentials (ip,a,j3) are replaced by the electron fluid velocity v. We also show 
that these two Lagrangian densities are in fact connected by a simple formula. Lastly, we 
summarize our work in Sec. V. 

II. CONSERVATION LAWS OF EULER-LAGRANGE EQUATIONS 

Standard Euler-Lagrange (EL) field equations are derived from the variational principle 
(JHJ) when the variations 5ip a are all independent. The variation 5C of the Lagrangian density 
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in Eq. (J2J) is, thus, explicitly written as 

dC 



5C 



E 



d 



dC 



d(d^ip a ) 



+ 



dx^ 



dip a dx^ 

where the exact space-time divergence is obtained by rearranging terms with 

dC 



A£l 



E 



(5) 



(6) 



did^ipa) 

The stationarity of the action functional associated with the variational principle Q yields 



= E / d 4 x5ip a 

„ JV 



dC 







dC 



d(d^ip a ) 



(7) 



iv y dip a dxv 

where the exact space-time divergence <9 M A M drops out under the assumption that the vari- 
ations 5ip a vanish on dT>. 

Following the standard rules of Calculus of Variations , the condition that Eq. (JJJ) holds 
for any variation 5ip a yields the Euler-Lagrange equation for the field-component ip a : 
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dx^ 



dC 



d(d^ip a ) 



d_ 

dt 



dC 



d(d t ip a ) 



+ V' 



dC 



d(Vijj a ) 



dC 

dipa 



(8) 



Since the Euler-Lagrange equations 
equation (jSJ) becomes 



hold for arbitrary variations 5ip a , the variational 



5C = 







E W« 



dC 



d(d^ip a ) 



(9) 



dx^ dx^ 

which is the Euler-Lagrange form of the Noether equation. 

The standard application of the Noether method involves the derivation of the energy- 
momentum conservation law. For this purpose, we consider arbitrary space-time translations 
[x v —>■ x v + 5x u ) generated by the space-time displacement 5x u . Under this space-time trans- 
formation, the variations of the field-components and the Lagrangian density are expressed 
as 

Slpa = - 5X V d v 1p a 

(10) 

8C = - d u (5x v C) + 5x v d' u C t 
where d' u C = {d v C)\$ a denotes the space-time derivative of C at constant field ip a . The 
energy-momentum conservation law is written in terms of the energy-momentum tensor 

dC 



E 



t 9{d^ip a ) 



d v 1p a - 



(11) 



as 



6^ v = d' v C. (12) 

Hence, when the Lagrangian density £ is independent of the space-time coordinate x x (i.e., 
d' x C = 0), Noether's theorem yields the energy-momentum conservation law 6^T\ = 0. 

When the background medium is time independent, the energy conservation law is written 
as 68 /dt + V • S = 0, where the energy density 8 and the energy-density flux S are defined 

as 
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6C 6i/j a 



6(6M dt 



— £ and S 



E 



6C 6i/j a 



6{V^ a ) dt 



(13) 



When the background medium is homogeneous, on the other hand, the momentum con- 
servation law is written as 611/ dt + V ■ T = 0, where the momentum density II and the 
momentum-stress tensor T are defined as 



n 



E 



dC 



6(6 t ip a ) 



and T 



£i - E 



6C 



6(Vip a 



Note that the energy conservation law is invariant under the transformation 

(s, s) - (e + v-c, s - 



(14) 



(15) 



where the vector C is an arbitrary function of space and time. The momentum conservation 
law, on the other hand, is invariant under the transformation 



in, t>~ (n + v.D,T-f) 



(16) 



where the tensor D is an arbitrary function of space and time. The invariance properties 
of the energy- momentum conservation law (T^, — > T^ u + d^K^, where K^ ieT = — so 
that cP = 0) can be used, for example, to ensure that the momentum-stress tensor T 
is symmetric. 



III. RELATIVISTIC COLD LASER-PLASMA INTERACTIONS 

The basic set of nonlinear equations describing the interaction of an intense laser pulse 
with a cold relativistic electron plasma [6| is given as 

6ti 

_ + V-nv = 0, (17) 
^ + v-Vp = -e (E + Ixb), (18) 
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VxB — 



V-E = 4vre (N - n) 
1 dE , v 




(19) 
(20) 



where p = Twy v is the relativistic electron kinetic momentum, with the relativistic factor 



while E = — Vy? — c~ 1 dtA and B = V X A are the electric and magnetic fields expressed in 
terms of the electromagnetic potentials (if, A). In Eq. ()19)) . the neutralizing fixed- ion density 
N is assumed to be nonuniform (VN ^ 0) but is taken to be time-independent (otherwise, 
the ion continuity equation would require ion flow). We note that, while the vector potential 
A includes both the laser field (Ao) and the field (A p ) induced by the plasma current 
response, the scalar potential <p represents only the field (y? p ) induced by the plasma charge- 
separation response. Furthermore, we note that the assumption of a cold relativistic electron 
plasma is justified by the large electron quiver velocities v g = — eAo/(m e c) generated by the 
intense laser field, which dominate over the random thermal motion of electrons. Although 
thermal effects can easily be incorporated within a Lagrangian formulation, they represent 
unecessary complications for our present purposes. 

We now show that this set of equations can be represented either in terms of an Euler- 
Lagrange variational formulation or an Euler-Poincare variational formulation. The Euler- 
Lagrange formulation presented below focuses its attention on the relativistic electron canon- 
ical momentum P = p — eA/c instead of the relativistic electron kinetic momentum p, which 
satisfies the relativistic canonical momentum equation 



obtained from Eq. (|18|). From Eq. (|22|). we note that the dynamical evolution of the mo- 
mentum vorticity O = V X P is expressed as dfl/dt = V X (v x O), which implies that 
the evolution of O requires that its initial value be nonzero (i.e., if the electron canonical 
momentum is initially curl- free, it remains so during its subsequent evolution). The Euler- 
Poincare formulation, on the other hand, considers variations in electron density n and fluid 
velocity v to be constrained by the conservation of mass, i.e., the variations 5n and <5v must 
satisfy the perturbed continuity equation 
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Both variational formulations yield identical energy-momentum conservation laws. 
The energy conservation law is, first, written in primitive form as 

1 _ , _ 1 



+ V- 



— (E + V^-E f|E| 2 - |B| 2 ) - e(n-N)(p + n 7 mc 2 

47T ' 8tt v ; 
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(24) 



Next, by removing the gauge vector C = E ip/An associated with the energy-gauge condition 
(fTHjl . we obtain 



d 
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n 7 mc 2 H (|E| 2 + IB I 
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(p e (n - N) + 
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V-E^ 

47T 



nv 7 mc + — E x B + ip — - — VxB - env 



All F ^4^ <9t 47T 

Lastly, making use of the Maxwell equations (|T9| and (J20|) and the continuity equation |T7j) 
we obtain the energy conservation law 

<9£ 







at 



v-s, 



where the energy density £ and energy-density flux S are 



E = n^mc 2 + — (JE| 2 + |B| 2 J , 
S = nv 7 mc 2 + — E X B. 

47T 

The momentum conservation is written in primitive form as 
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We now remove the gauge tensor D = E A/47rc associated with the momentum-gauge con- 
dition (fTT)j) to obtain 
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Next, making use of the Maxwell equations (|19|) and (|20|). we obtain 



ep ViV 
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' dt 
+ V- 



n m7v + 
m7 vv 
which finally becomes 



E x B 
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where the momentum density II and the symmetric momentum-stress tensor T are 

E x B 
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T 



n m'-fv + 



nm'-f vv 
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iTT^ ' 8tt 



El 
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(29) 

(30) 
(31) 



Hence, Eq. ()29|) implies that the total momentum is conserved in the direction where the 
electric field vanishes. Note that, while the momentum-stress tensor in the primitive form 
of the momentum conservation law (}2"%|) is not symmetric, the final form of the momentum- 
stress tensor (}3lf is symmetric. 



A. Euler-Lagrange Formulation 

The Euler-Lagrange formulation of the relativistic cold laser-plasma equations (fTTj) - (l2"Ul) is 

based on the use of Clebsch (or Euler) potentials p| . According to the Clebsch representation 
associated with the general vorticity condition (O = V X P ^ 0), the relativistic electron 
canonical momentum is expressed as 

P = m7 v - - A = + a V/3, (32) 

where the Clebsch potentials ip, a, (3 are functions of space and time. As will be shown 
below, the potential ip may be viewed as a Lagrange multiplier associated with the mass 
conservation law (e.g., the continuity equation), while the potentials a and (3 ensure that 
V X P = V« X V/3 7^ 0. Furthermore, the Lin constraints jj| on the Clebsch potentials a 
and P requires that the potentials a and (3 be carried by the fluid flow in the sense that 



^ + v-V« = = f t + vVft (33) 

i.e., the potentials a and (3 are known as Lagrangian- marker coordinates 0,1^. Note that 
under a gauge transformation ((p, A) — » (jp — c~ x d t \i A + Vx), the Clebsch potential ip 



obeys the gauge transformation ip — > ip — (e/c) x while the Lagrangian-marker coordinates 
a and (3 are invariant. Lastly, a more symmetric expression for the Clebsch representation 
of the relativistic canonical momentum can also be obtained by writing Eq. (|32j) as 



Vip' + - (orV/3 - /3Va) 



(34) 



which involves a redefinition of the Clebsch potential ip' = if) + ~ a/3. 

The Lagrangian density in the Euler-Lagrange formulation of the relativistic cold laser- 
plasma equations {HJ), (HH-flU, and (j22J) is 
1 
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where the relativistic factor is 
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(35) 



(36) 



where the relativistic canonical momentum is expressed in terms of Clebsch potentials as 
given in Eq. (}3"2j) . Here, the variational fields (ip, A; n, if), a, f3) are all to be varied indepen- 
dently in the variational principle 



5 I Cel(I a x 



0. 



(37) 



The variation of the Lagrangian density (|33j) can be expressed as 
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(38) 



where we have rearranged terms in order to isolate variations in the dynamical variables 
(<p, A; n, ip, a, (3) and, thus, extracted the Noether space-time divergence. 

As a result of the variational principle (J3"T|) . the relativistic cold plasma-laser equations 
(fTTj) . (p"9|) - (j20|) . and (J33|) are easily recovered. The last variational equation, which corre- 
sponds to the Euler-Lagrange equation dCEi/dn = 0, is the energy equation 

,2 



dip dp 



e cp 



7 mc 



(39) 



We recover the relativistic electron canonical momentum equation (}22j) by first taking the 
gradient of Eq. and rearranging terms to obtain 

d t {Vip + aV/3) = S7 (&p - ^mc 2 ^ + vx (Va X V/3) . 

By substituting the Clebsch representation we easily obtain Eq. Since these 

variational equations hold for arbitrary variations in (ip, A; n, ip, a, (3), the expression for the 
variation (J38)) of the Lagrangian density can now be written as the Euler-Lagrange form of 
the Noether equation 
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(40) 



which can now be used to derive the conservation laws of energy and momentum for the 
relativistic cold laser-plasma equations. 

1. Energy conservation law for the relativistic cold laser-plasma equations 



The conservation law of energy is derived from the Noether equation (J4*U|) by considering 
time translations t — > t + 5t, such that Sip = — 5t d t ip, 5(3 = — 5t d t (3, 5p> = — 5t d t <p, 

5A = -5td t A = cSt(E + Vy>), 

and SCel = — 5td t £EL (where we use the fact that the fixed-ion density iV is time- 
independent). Substituting these expressions into Eq. (|4T)j) . we first obtain 
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V <9t dt / 47r dt 4n 



which, by making use of Eq. (J39|) . becomes the primitive form of the energy conservation 
law and, thus, yields the energy conservation law ([25)1 . 



Momentum conservation law for the relativistic cold laser-plasma equations 



The conservation law of momentum is derived from the Noether equation (J40j) by con- 
sidering space translations x — > x + 5x, such that Sip = — 5x • Vtp, 5(3 = — <5x ■ V/3, 

1 0A> 
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and SCel = — V ■ (<5x Cel) + Sx-VCel (where V Cel = — eip ViV for a nonuniform 
fixed-ion density). Substituting these expressions into Eq. (}4Tjj) . we first obtain 

1 
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Next, substituting the Clebsch representation (|H2|) for the relativistic canonical momen- 
tum, we obtain the primitive form (|28|) of the momentum conservation law and, thus, the 
momentum conservation law ()29|) follows. 



B. Euler-Poincare Formulation 



The Euler-Poincare variational formulation of the relativistic cold laser-plasma equations 
flT7j) - (j2T?|) is expressed in terms of the Lagrangian density 
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8tt v ' ' c 

where 7~ x = (1 — |v| 2 /c 2 )2. The Eulerian variations of the electron density n and fluid 
velocity v, which are subject to the mass conservation constraint (J2HJ), are expressed in 
terms of the virtual fluid displacement 5$, as 
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v • V )5£ - 5£ • Vv. 



(42) 



Relativistically covariant expressions for the Eulerian variations for the proper density iV = 
7 _1 n and the fluid four- velocity = (70,7V) are found in Ref. here, we prefer the 
Eulerian variations of n and v for the sake of simplicity. 

The variation of the Euler-Poincare Lagrangian density can be expressed as 
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where P = rwyv — (e/c) A denotes the relativistic canonical momentum and we have rear- 
ranged terms in order to isolate variations in the dynamical variables (£, tp, A) and, thus, 
extracted the Noether space-time divergence. 
As a result of the variational principle 

S f C EP d A x = 0, (44) 



the relativistic cold laser-plasma equations (fTT)|) - (|2T)|) and (j22J) are easily recovered (under 
arbitrary variations Sip, SA, and 5£, respectively) with the continuity equation (jl7j) repre- 
senting a constraint equation. Since these variational equations hold for arbitrary variations 
in (£,<£>, A), the expression for the variation (j4Hj) of the Lagrangian density can now be 
written as the Euler-Poincare form of the Noether equation 



SC EP = - -^-SA-E 

at \ A-kc 



nS£ (7 1 mc 2 — ap + -v-A 



+ nvP-££ - —ISuE + SAxB 

4tt 



(45) 



which can now be used to derive the conservation laws of energy and momentum for the 
relativistic cold plasma-laser equations. 

To derive the energy-momentum conservation laws within the Euler-Poincare formulation, 
the virtual fluid displacement 5£ is expressed in terms of the space and time translations #x 
and St as 

S£ = 5x - vSt, (46) 

where v is the fluid velocity. Note that the relativistic approach used here is based on con- 
strained variations of the laboratory-frame electron fluid density and electron fluid velocity; 
an alternative approach is based on constrained variations of rest-frame fluid quantities. 

1. Energy conservation law for the relativistic cold laser-plasma equations 

The conservation law of energy is derived from the Noether equation (}45|) by considering 
time translations t — > t + St, such that S£ = — v St, Sip = — St d t (p, 

SA = -Std t A = cSt(E + V<p), 
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and SCep = — Std t CEP (where we use the fact that the fixed-ion density N is time- 
independent). Substituting these expressions into Eq. (J4~5j) . we first obtain 
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Next, by substituting the expression 
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we recover the primitive form (}24|) of the energy conservation law, from which the energy 
conservation law follows. 



2. Momentum conservation law for the relativistic cold laser-plasma equations 



The conservation law of momentum is derived from the Noether equation (j4*3j) by con- 
sidering space translations x — > x + <5x, such that 5£ = 5x, 



6(p = - <5x • Vy? = 5x • E H — , 

V c at } 



5A 



5x-VA = fexB - V(A-<5x) 



and 5Cep = — V ■ (5x £_ep) + 5x-V£ep (where VCep = — eip ViV for a nonuniform 
fixed-ion density). Substituting these expressions into Eq. (|43|) . we first obtain 
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Next, by substituting the expression P = 777/yv — (e/c) A for the relativistic canonical mo- 
mentum, we recover the primitive form (|2K|) of the momentum conservation law, from which 
the momentum conservation law ()29|) follows. 



C. Connection between the Lagrangian Densities 

We have thus far shown that the relativistic cold laser-plasma equations (fTTjl - pUJ) can 
be represented in terms of two different variational principles based on either the Euler- 
Lagrange Lagrangian density (J33j) or the Euler-Poincare Lagrangian density ([41)1 . We now 
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show the connection between these two Lagrangians. We begin by evaluating the difference: 
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Next, by making use of Eq. ()47|) with the Clebsch representation (J3"4^ for the canonical 
momentum P, the difference (li'81) becomes 
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Here, the first term can be written as an exact space-time divergence plus a term involving 
the continuity equation and, thus, the difference Cel — Cep is expressed as 
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(49) 



Lastly, since two Lagrangian densities that differ by an exact space-time divergence are 
considered equivalent (i.e., they generate identical dynamical equations), we may eliminate 
the exact space-time divergence from the right side of Eq. (j4Uj) to obtain the final expression 



£el = £ep + *l>' 



dn 
~di 



+ V-(nv) - 



n 



v • Va 



(50) 

for the Euler-Lagrange Lagrangian density £el given in terms of the Euler-Poincare La- 
grangian density £ep, the mass-conservation constraint (fTTjl (with Lagrange multiplier ip'), 
and the Lin constraints (|33j) for the Lagrangian-marker coordinates a and (3 (with Lagrange 
multipliers n/3/2 and —na/2, respectively). 



IV. SUMMARY 



We have presented two variational formulations for the relativistic cold laser-plasma equa- 
tions, which depend on whether all the dynamical fields are varied independently or not. 
We note that variational formulations in the physics literature are either presented in the 
Euler-Lagrange form or the Euler-Poincare form (generally not both) and that the work 
presented here shows the explicit connection ()50|) between the two variational formulations 
for an important set of plasma-fluid equations. 
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In deriving the energy-momentum conservation laws for the relativistic cold laser-plasma 
equations, we showed that, although the form of the Noether equation depends on the 
variational formulation used, i.e., either the Euler-Lagrange form (J4U)) or the Euler-Poincare 
form (J45j) . the same conservation equations (J25j) and (|29|) are obtained. 
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